Coherent pulsed excitation of degenerate multistate systems: Exact analytic solutions 
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We show that the solution of a multistate system composed of N degenerate lower (ground) 
states and one upper (excited) state can be reduced by using the Morris-Shore transformation to 
the solution of a two-state system involving only the excited state and a (bright) superposition of 
ground states. In addition, there are A'^ — 1 dark states composed of ground states. We use this 
decomposition to derive analytical solutions for degenerate extensions of the most popular exactly 
soluble models: the resonance solution, the Rabi, Landau-Zener, Rosen-Zener, Allen-Eberly and 
Demkov-Kunike models. We suggest various applications of the multistate solutions, for example, 
as tools for creating multistate coherent superpositions by generalized resonant 7r-pulses. We show 
that such generalized 7r-pulses can occur even when the upper state is far off resonance, at specific 
detunings, which makes it possible to operate in the degenerate ground-state manifold without 
populating the (possibly lossy) upper state, even transiently. 

PACS numbers: 32.80.Bx, 32.80.Qk, 33.80.Be, 32.80.-t, 33.80.-b 



I. INTRODUCTION 

The problem of a two-state quantum system driven 
by a time-dependent pulsed external field plays a central 
role in quantum physics First of all, this problem is 
interesting by itself both physically and mathematically: 
physically, because the two-state system is the simplest 
nontrivial system with discrete energy states in quan- 
tum mechanics; mathematically, because the Schrodinger 
equation for two states poses interesting mathematical 
challenges some of which are exactly soluble. Further- 
more, already in the two-state case, important nonclassi- 
cal phenomena occur, for instance, the famous Rabi oscil- 
lations, which often serve as a test for quantum behavior, 
and also provide a powerful tool for coherent control of 
quantum dynamics, e.g. by tt pulses. Finally, in almost 
all cases (except for a few exactly soluble), the behavior 
of a multistate quantum system can only be understood 
by reduction to one or more effective two-state systems, 
e.g., by adiabatic elimination of weakly coupled states or 
by using some intrinsic symmetries. 

Besides the well-known solution for exact resonance, 
there exist several exactly soluble two-state models, the 
most widely used being the Rabi 0, Landau-Zener 3|, 
Rosen-Zener Allen-Eberly [sl], Bambini-Berman ^, 
Demkov-Kunike 0, CarroU-Hioe [S], Demkov Q and 
Nikitin [T^j models. All these models provide the transi- 
tion probability between two nondegenerate states. 

In the present paper, we present the extensions of these 
exactly soluble models to the case when one of the states 
is replaced by TV degenerate states, as displayed in Fig . [1] 
By using the Morris-Shore (MS) transformation [ll| we 
show that the (A^-f l)-state problem can be reduced to an 
effective two-state problem involving a bright state and 
the upper, nondegenerate state. If known, the propaga- 
tor for this subsystem can be used to find the solution for 
the full (iV-l- l)-state system. Such analytic solutions can 



be very useful in designing general unitary transforma- 
tions within the iV-state degenerate manifold, which can 
be viewed as a qunit for quantum information process- 
ing We point out that the same system for = 3 
has been considered by Unanyan et al [l^ and by Kis 
and Stenholm lj| for general N, who have derived the 
adiabatic solution for pulses generally delayed in time; 
these schemes extend the well-known technique of stim- 
ulated Raman adiabatic passage (STIRAP) (see for 
reviews) . Here we derive several exact analytic solutions 
for pulses coincident in time. This work can therefore be 
considered as an extension to arbitrary N of an earlier 
paper jT3|, which treated the case N = 2. 

This paper is organised as follows. In Sec. [IT] we de- 
scribe the system and define the problem. In Sec. IIIII 
we introduce the MS basis and derive the {N + l)-state 
propagator in terms of the (presumably known) two-state 
propagator. In Sec. |TV] we use this solution to iden- 
tify various interesting types of population evolutions. 
In Sec. |V] we use the analytic solutions for exact res- 
onance and the Rosen-Zener model to propose several 
applications, for example, creation of maximally coher- 
ent superpositions and qunit rotation. In Sec. IVII we 
discuss some aspects of the multistate Landau-Zener and 
Demkov-Kunike models. Finally, Sec. IVIII provides a 
summary of the results. 



II. DEFINITION OF THE PROBLEM 

A. System Hamiltonian 

We consider an [N + l)-state system with N degener- 
ate lower (ground) states l^/^n) {n = l,2,...,iV) and one 
upper (excited) state \ij}N+i), as depicted in Fig. [TJ 
The N lower states are coupled via the upper state 
with pulsed interactions, each pair of which are on two- 
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FIG. 1: Top: The system studied in this paper. A*' degen- 
erate (in RWA sense) states It/ji) , |i/'2) , . . . , IV'iv) ^^re coupled 
simultaneously to an upper state |?/>jv+i), possibly off single- 
photon resonance by a detuning A(t), with Rabi frequencies 
f^n(i) (n — 1,2, . . . , N). Bottom: The same system in the 
Morris-Shore basis. There are A*' — 1 uncoupled dark states 
\fi) , \^P2) , ■ ■ ■ , j<^jv-i), and a pair of coupled states, a bright 
state Iv^at) and the upper state \ipN+i), with the same detun- 
ing A{t) as in the original basis and a coupling given by the 
rms Rabi frequency fl{t), Eq. 



photon resonance (Fig. [T]). The upper state \ipN+i) 
may be off single-photon resonance by some detuning 
A{t) that, however, must be the same for all fields. 
In the usual rotating-wave approximation (RWA) the 
Schrodinger equation of the system reads [l| 



ih-Cit) = H(i)C(t), 



(1) 



where the elements of the {N + l)-diniensional vector 
C{t) are the probability amplitudes of the states and the 
Hamiltonian is given by 







ill (t) 

^2 {t) 



nM{t) 

ilNit) 2A{t). 



(2) 



For the sake of simplicity the Rabi frequences of the cou- 
plings between the ground states and the excited state 
V,i {t) , {t) are assumed real and positive as the 

populations do not depend on their signs. The phases of 
the couplings can easily be incorporated in the descrip- 
tion and they can be used to control the inner phases 
of the created superposition states. Furthermore, the 



Rabi frequencies are assumed to be pulse-shaped func- 
tions with the same time dependence f{t), but possibly 
with different magnitudes. 



(t) = Xnf{t) (n : 

and hence different pulse areas. 



1,2, 



(3) 



An — 



Vln{t)dt = Xn 



f{t)dt (n=l,2,...,iV). 

J 

(4) 



B. Physical implementations 

The linkage pattern described by the Hamiltonian ([2]) 
can be implemented experimentally in laser excitation of 
atoms or molecules. For example, the iV = 3 case is 
readily implemented in the J = 1 <-> J = system cou- 
pled by three laser fields with right circular, left circular 
and linear polarizations, as shown in Fig. [2] (left). These 
coupling fields can be produced from the same laser by 
standard optical tools (beam splitters, polarizers, etc.), 
which greatly facilitates implementation. Moreover, the 
use of pulses derived from the same laser ensures auto- 
matically the two-photon resonance conditions and the 
condition ([3]) for the same temporal profile of all pulses. 

The cases of iV = 4 — 6 can be realized by adding 
an additional J = 1 level to the coupling scheme and 
appropriately polarized laser pulses, as shown in the right 
frame of Fig. [2l 
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FIG. 2: Examples of physical implementations of the link- 
age pattern of N degenerate ground states coupled via one 
upper state, considered in the present paper. Left: A'^ = 3 
degenerate states. Right: A'' = 4 degenerate (in the RWA 
sense) states (dashed arrows indicate two additional possible 
linkages) . 



III. GENERAL SOLUTION 

A. Morris-Shore (dark-bright) basis 

The Hamiltonian ([2|) has — 1 zero eigenvalues and 
two nonzero ones. 



A±(t) 



A„ = (n = 1 
1 r 



,N-l) 



(5a) 
(5b) 



3 



where 



N 



\ n=l 



(6) 



is the root-mean-square (rms) Rabi frequency, where 



N 



\ n=l 



(7) 



The set of orthonormaUzed eigenstates \(pn) {n = 
l,2,...,iV— 1) corresponding to the zero eigenvalues can 
be chosen as 



\fN-l) 

where 



1 T 
-TT [X2,-Xi,0,0,- • • ,0] , 

^2 



(8a) 



^ XiX3,X2X3,-^|,0,--- ,0]^, (8b) 



X2X3 

y V [xiX4, X2X4, X3X4, -X^, 0, • • • , 

A3A4 



1 



Xn-iXn 



[xiXN, X2XN, ■■■ , -x^„i, o](i4) 



Xn — 



k=l 



Xl (n-2,3,. 



(9) 



These eigenstates are dark states, i.e. they do not involve 
the excited state \iI>n+i) and, as we shall see, are uncou- 
pled from IV-'Tv+i). All dark states are time-independent. 
We emphasize that the choice ^ of dark states is not 
unique because any superposition of dark states is a dark 
state too; hence their choice is a matter of convenience. 

The Hilbert space is decomposed into two subspaces: 
an [N — l)-dimensional dark subspace comprising the 
dark states ([5]) and a two-dimensional subspace orthog- 
onal to the dark subspace. It is convenient to use the 
Morris-Shore (MS) basis which, in addition to the 
dark states, includes the excited state = Wn+i) 

and a bright ground state \^n)- The latter does not have 
a component of the excited state and is orthogonal to the 
dark states; these conditions determine it completely (up 
to an unimportant global phase). 



1 5 
\vn) = — [xi,X2, ■ • • ,Xw,0] 



(10) 



We point out that the Morris-Shore basis is not the adi- 
abatic basis because only the dark states are eigenstates 
of the Hamiltonian, but I^at) and |</JAr+i) are not. 

In the new, still stationary basis {\^n)}n= 
the Schrodinger equation reads 



ihjB{t) = H(i)B(i), 



(11) 



where the original amplitudes C{t) are connected to the 
MS amplitudes B(i) by the time-independent unitary 
matrix W composed by the basis vectors \(pn), 



\N=[\vi),\V2) 



according to 



C{t) = \NB{t). 



(12) 



(13) 



The transformed Hamiltonian reads H(t) = Wt|-|(t)W, or 
explicitly. 



-0 


• 


• 








• 


• 








• 


• 








• 


• 


n{t) 


.0 


• 


• n{t) 


2A{t) 



(14) 



We point out that the Hamiltonian of Eq. ((21) is a 
special case of the most general Hamiltonian for which 
the MS transformation [ll| applies and which includes N 
degenerate lower states and M degenerate upper states. 
Hamiltonians of the same type as ^ and related trans- 
formations leading to Eq. (HH), have appeared in the 
literature also before the paper by Morris and Shore [ll| , 
mostly in simplified versions of constant and equal inter- 
actions (see e.g. [l^ and references therein). 



B. Solution in the Morris-Shore basis 

As evident from the first — 1 zero rows of H the dark 
states are decoupled from states \(Pn) and \ipN+i) and the 
dark-state amplitudes remain unchanged, Bn{t) — const 
{n = 1,2, . . . ,N - 1). Thus the (A^ + l)-state problem 
reduces to a two-state one involving \ipN) and 



dt 
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'0 


n ' 




Bn 


Bn+1 


~ 2 


n 


2A 




Bn+1 



(15) 



The propagator for this two-state system, defined by 



Bn (+00) 
Bn+1 (+00) 



Bn (-00) 
Bn+1 (-00) 



(16) 



is unitary and can be expressed in terms of the Cayley- 
Klein parameters as 



U 



(2) 
MS 



a 

-b* 



b 

a* 



(17) 



with \b\^ = 1 — |a|^. Then the transition matrix for the 
(A^ -I- l)-state system in the MS basis reads 



'MS 



1 
1 





.0 





a 

-b* 




b 

a* 



(18) 
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C. The solution in the original basis 



or explicitly, 



We can find the transition matrix in the original, dia- 
batic basis by using the transformation 

U(^+iHoo,-oo) = WUi^/'^(oo,-c5o)Wt, (19) 



(a-1) 
1 + (a - 1 
(a-1) 



(a 



1) XiXiV 
X 



XlX2 

X2 
X2X3 



(a-1) 



X2XiV 

5* X2 

X 



1) X1X3 

^ x^ 



(a 
(a 

1 + (a - 1) 



X2X3 

x:^ 
x^ 



(a-1) 



X3XiV 
5*X3 

X 



a — 


1) 


XlXJV 


5Xi 






x' 


X 


a — 


1) 


X2XN 








X' 


X 


a — 


1) 


X3XJV 


&X3 






x^ 


X 



1 + (a - 1) ^ 6^ 

X 



(20) 



The ith column of this matrix provides the probability 
amplitudes for initial conditions 



Q(-(^) = 1, 
C„(-oo) - (n^i). 



(21a) 
(21b) 



The initial state 1-04) can be one of the degenerate states 
or the upper state. This general unitary matrix and com- 
binations of such matrices can be used to design tech- 
niques for general or special qunit rotations. 

As evident from Eq. for finding the populations 
for the initial condition pT|) it is sufficient to know only 

the parameter a 

\a!\ For the sake of simplicity, in the present paper 

we are interested only in cases when the system starts 
in a single state and below we shall concentrate on the 
values of the parameter a. In the more general case when 
the system starts in a coherent superposition of states, 
Eq. (pn|) can be used again to derive the solution; then 
the other Cayley-Klein parameter h is also needed. 



C/|^g(oo, — oo) because \h\^ — 1 ■ 



(PO]) that the populations in the end of the evolution are 



P. 
P.. 



1 + (a - 1) 



la- 1 



2 2 
|2 Xn 

4 

X 



Pn+1 = (1-H') ^. 



(22a) 

(n^i,7V+l), (22b) 
(22c) 



Therefore the ratio of the populations of any two degen- 
erate states, different from the initial state \tpi), reads 



Pm _ Xm 
Pn ^ Xl 



2 (m,n ^ i,7V+ 1). 



(23) 



Hence these population ratios do not depend on the inter- 
action details but only on the ratios of the corresponding 
peak Rabi frequencies. 

For equal Rabi frequencies. 



IV. TYPES OF POPULATION DISTRIBUTION 

We identify two types of initial conditions: when the 
system starts in one of the degenerate states or in the 
excited state |?/'Ar_|_i), which we shall consider separately. 



A. System initially in a ground state 

1. General case 



Xl = X2 = ■ ■ ■ = XN, 

Eqs. ([22]) reduce to 

2 



P = 



Pn = 



1 



a-1 



P, 



N+l 



\a~ 1, 
N 



N 

2 



(24) 



(25a) 

n^i,N+l), (25b) 
(25c) 



When the system is initially in the ground state 
Eq. (|2ip . we find from the ith column of the propagator 



Thus the populations of all ground states except the ini- 
tial state \ipi) are equal. 
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2. Special values of a 



if 



Several values of the propagator parameter a are espe- 
cially interesting. 

For a = 0, which indicates complete population trans- 
fer (CPT) in the MS two-state system, Eq. ^ gives 



P. 

Pn 

Pn+1 



An ^2 
,.2 



(26a) 
(26b) 
(26c) 



For a — 1, which corresponds to complete population 
return (CPR) in the MS two-state system, we obtain 



P = 1, 

Pn = {n^i,N + l), 
Pn+1 = 0. 



(27a) 
(27b) 
(27c) 



For a — —1, which again corresponds to CPR in the 
MS two-state system, but with a sign flip in the ampli- 
tude, we find 



P. = 1-2^ 



(28a) 



Pn = {n^t,N+l), (28b) 



N+l 



= 0. 



(28c) 



It is important to note that although both cases a = 1 
and a = — 1 lead to CPR in the MS two-state system, 
they produce very different population distributions in 
the full {N + l)-state system. The case a = 1 leads to a 
trivial result (CPR in the full system), whereas the case 
a = — 1 is very interesting because it leads to a popula- 
tion redistribution amongst the ground states with zero 
population in the upper state; hence this case deserves a 
special attention. 



3. The case a — —1 

The case of a = —1 is particularly important because 
it allows to create a coherent superposition of all ground 
states, with no population in the upper state. 

All ground-state populations in this superposition will 
be equal. 



Pl=P2 = 

Pn+1 = 0. 



P 



N 



1 

iV' 



(29a) 
(29b) 



where 



Xn = 

Xo = 



(%/]V±l)xo, 



(30a) 
(30b) 

(30c) 



This result does not depend on other interaction details 
(pulse shape, pulse area, detuning) as long as a = — 1. 
For example, for = 4 degenerate states, equal popula- 
tions are obtained when Xi — Xn or Xi = 3x„. We shall 
discuss later how the condition a — —I can be obtained 
for several analytically soluble models. 

Another important particular case is when the initial- 
state population Pi vanishes in the end. This occurs for 



X^ 



(31) 



For example, an equal superposition of all lower sublevels 
except IV'i), 



P, = Pn+1 = 

p.. ^ ' 



N- 1 



is created for 



Xi 

Xn 



= XoVN^, 
= Xo {n^ i), 



where 



Xo 



^2(^-1) 



(32a) 
(32b) 



(33a) 
(33b) 



(33c) 



B. System initially in the upper state 

If the system is initially in the excited state j'^jv+i) 
the end of the evolution the populations will be 



at 



Pn=[l-\a\^)^ ("-1,2,, 



Pna 



a 



,N), (34a) 
(34b) 



For a = ±1 at the end of the evolution the system un- 
dergoes CPR, as in the MS two-state system. For a = Q 
(CPT in the MS two-state system) the whole popula- 
tion will be in the ground states leaving the excited state 
empty, Pn+i = 0. If all the couplings are equal, Eq. 
I, the ground states will have equal populations. 



Pn 



N ^ 



,N). 



(35) 
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C. Discussion 

In this section we discussed some general features of 
the population redistribution in the (Af + l)-state system. 
There are three particularly interesing results. 

First, the ratios of the populations of the degenerate 
states (except the one populated initially) depend only on 
the ratios of the corresponding Rabi frequencies; hence 
they can be controlled by changing the corresponding 
laser intensities alone. The populations values, though, 
depend on the other interaction details. Moreover, it can 
easily be seen that the relative phases of the degenerate 
states can be controlled by the relative laser phases. 

Second, it is possible to create an equal superposition of 
all ground states, with zero population in the upper state. 
This is possible when the system starts in a ground state: 
then condition (jSO]) is required, along with the CPR con- 
dition a = — 1 . Alternatively, an equal superposition can 
be created when the system starts in the upper state: 
then condition is required, along with the CPT con- 
dition a = 0. Equal superpositions are important in some 
applications because they are states with maximal coher- 
ence (since the population inversions vanish). 

Third, it is possible, starting from a ground state, to 
create a superposition of all other ground states, whereas 
the initial ground state and the excited state are left un- 
populated. This requires a = —1 and condition (j3ip . 
This case has interesting physical implications, which will 
be discussed in the next section. 



APPLICATIONS TO EXACTLY SOLUBLE 
MODELS 

A. Multistats analytical solutions 



The values 

r(2) 



of the propagator parameter a — 
U^'g{oo,—oo) for the most popular analytically ex- 
actly soluble models are listed in Table HI Equation ((20)) . 
supplied with these values, provides several exact multi- 
state analytical solutions, which generalize the respective 
two-state solutions. 

Among these solutions, the resonance case is the sim- 
plest and most important one, which will receive a special 
attention below. It will be followed by a detailed discus- 
sion of the Rosen-Zener (RZ) model, which can be seen as 
an extension of the resonance solution to nonzero detun- 
ing for a special pulse shape (hyperbolic secant). Both 
the resonance and the RZ model allow for the param- 
eter a to obtain the important values 0, ±1. The Rabi 
model can also be used to illustrate the interesting cases 
of population distribution associated with these values of 
a but its rectangular pulse shape is less attractive (and 
also less realistic) than the beautiful sech-shape of the 
pulse in the RZ model. 

The Landau-Zener (LZ) and AUen-Eberly (AE) mod- 
els are of level-crossing type, i.e. the detuning crosses 



Model 



Resonance 
a = cos 



Rabi 

m = xm^T), A{t) = Ao 

a = COS (t^x^ + - i ^ 



Landau-Zener 

n(t) = X, A(t) = ct 

a = exp [-7rxV4C] 



Rosen-Zener 

n{t) = xsech(t/r), A(t) = Ao 



AUen-Eberly 

n{t) = xsech(t/r), A{t) = Btanh(t/r) 



cosh (ttP) 



Demkov-Kunike 

n{t) = xsech(t/r), A(t) = Ao + Btanh(t/T) 

_ ^(i+^(^ + /3))^(i + ^(^-/3)) 



a2 - /32 + iS 



TABLE L Values of the Cayley-Klein parameter a 



U 



(2) 



J^^g(oo, — oo)j for several exactly soluble models. Here 

T{z) is the Gamma function and a — \xT, P = \BT, 
5 = iAoT, are scaled dimensionless parameters, which are 
assumed positive without loss of generality. 



resonance, A(0) = 0. For these models in the adiabatic 
limit the transition probability approaches unity, that is 
a ^ 0. The parameter a is always nonnegative, i.e. the 
most interesting value in the present context, a = — 1, is 
unreachable. Nevertheless, because of the popularity and 
the importance of the LZ model, and because the present 
multistate LZ solution supplements other multistate LZ 
solutions, we discuss this solution in detail in Sec. IVII 

The Demkov-Kunike (DK) model is a very versatile 
model, which combines and generalizes the RZ and AE 
models. Indeed, as seen in Table U the DK model re- 
duces to the RZ model for B = Q and to the AE model 
for Aq = 0. For the DK model, the parameter a can 
be equal to the most interesting value of —1 only when 
5 = 0, i.e. only in the RZ limit. Therefore, we shall 
only consider the RZ model below, and leave the AE and 
DK models to readers interested in other aspects of the 
analytic multistate solutions presented here. 



B. Exact resonance 



In the case of exact resonance, 
A = 0, 



(36) 
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the elements of the evolution matrix for the MS two-state where A is the rms pulse area defined as 
system for any pulse shape of ^l{t) are 



cos 



2 ' 



K ■ ■ ^ 

= — zsm — , 

2 ' 



(37a) 
(37b) 



A= / n{t')dt'. 

'J —OO 

In the important case of iV = 3 we have 



(38) 



I (4) 



-2^ 

X1X2 
x^ 



,2 

,2 _ 



-220X2 . 



2 1 



-2f*sm 4 

.22i«^ sin^ iA 

x^ 4 
_iX2 1^ 

X 2 



.220X3 sij^: 

X'' 

O X2X3 

x'-^„ 



1 - 2 



X3 



sm 



-i^ sm 
X 2 

-i^ sin 

X 2 

-i20 sin iyl 
X 2 

cos iv4 



(39) 



We have a = 0, ±1 for the following pulse areas, 



Multistats Rosen-Zener model 



a = 
a= 1 
a = -1 



A = (2? + 1) TT, (40a) 
A = 4:Itt, (40b) 
A = 2(2; + l)7r. (40c) 



where Z = 0, 1, 2, .... 

The pulse areas for the three important cases discussed 
in Sec. llVCI are easily calculated. 

An equal superposition of all N ground states is cre- 
ated when starting from the excited state and all indi- 
vidual pulse areas are equal to (see Sec. IIVB|) 

^^^^{21+^ (n = l,2,...,iV), (41) 



N 



where I = 0,1, 2, .... 

An equal superposition of all N ground states is cre- 
ated also when starting from one ground state and 
the pulse areas are [see Eq. (|5n|) ] 



A, 



An — 



'N ±1 
12 '^ {2l + l)7T, 



N 



(42a) 



N±VN 



{21 + 1)tt (n^i), (42b) 



where / = 0, 1, 2, ... 

The other interesting case when the system starts in 
one ground state and ends up in an equal superposi- 
tion of all other ground states is realised for pulse areas 
[see Eq. 



(43a) 



A, = \/2(2/-f l)7r, 

An = 



TV- 1 



(2? + l)7r (n^i), (43b) 



where Z = 0, 1, 2, 



Equation (|20p and the value of the parameter a in Ta- 
ble U represent the multistate RZ solution in the degen- 
erate two-level system. It is easy to show that 



= 1 - 



cosh^ (^ttAoT) ' 



(44) 



where we have used the reflection formula r(i + z)r(i — 
z) — tt/costtz T8|. Hence in this model |a| = 1 for 
a — = I {I = 0, 1,2, ...). The phase of a, however, 
depends on the detuning Aq fl^; we use this to an ad- 
vantage to select values of Aq for which a = —1. For 
a = I we find 



\j2l + l-tAoT 



k=0 



2l + l + iAaT 



(45) 



where the recurrence relation r(z + 1) = zT{z) ^T§\ has 
been used. Thus, the equation a — — I reduces to an alge- 
braic equation for Aq, which has I real solutions [17]. The 
first few values of x a-nd Aq for which a = — 1 are shown 
in Table HIl As the table shows, Aq = is a solution for 
odd a = ^xT but not for even a, in agreement with the 
conclusions in Sec. IV Bl Moreover, the a = — 1 solutions 
do not depend on the number of degenerate states N. 

In the present context the RZ model is interesting for 
it shows that one can create superpositions within the 
ground-state manifold even when the excited state is off 
resonance by a considerable detuning (Aq ^ ^/T), for 
which the transition probability in the MS two-state sys- 
tem is virtually zero, i.e. jaj w 1. This fact allows us, for 
specific detunings, to essentially contain the transient dy- 
namics within the ground states; in contrast, in the reso- 
nance case the excited state can get significant transient 
population, PN+i{t) = sin^ 5^(^)1 although it vanishes 
in the end. 

Figure [3] displays the populations against the detun- 
ing Aq for a hyperbolic-secant pulse with xT — 18 for 
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AoT 


2 















4 


±1.732 












6 





±4.796 










8 


±1.113 


±9.207 










10 





±2.756 


±14.913 








12 


±0.943 


±4.936 


±21.903 








14 





±2.243 


±7.595 


±30.171 






16 


±0.855 


±3.916 


±10.708 


±39.715 






18 





±1.988 


±5.907 


±14.265 


±50.534 




20 


±0.799 


±3.418 


±8.195 


±18.260 


±62.627 




22 





±1.830 


±5.098 


±10.766 


±22.687 


±75.993 


24 


±0.759 


±3.113 


±7.006 


±13.613 


±27.545 


±90.634 


26 





±1.719 


±4.606 


±9.130 


±16.729 


±32.833 






±106.549 










28 


±0.728 


±2.901 


±6.289 


±11.461 


±20.113 


±38.548 






±123.736 










30 





±1.636 


±4.268 


±8.150 


±13.994 


±23.760 






±44.690 


±142.198 









TABLE II: Some approximate solutions of the equation 
a(Ao) = —1 for the RZ model, where a is given in Table 
|I1 for various even integer values of xT. 



couplings chosen to satisfy Eqs. (|30p (upper frame) and 
([55)1 (lower frame). In both cases we have |a| = 1 [see 
Eq. (|44p ]. which leaves the excited state unpopulated in 
the end. For several special values of the detuning Aq, as 
predicted in Table HIl we have a = — 1. For these values, 
an equal superposition of all degenerate states including 
the initially populated state IV"!) is created in the upper 
frame, and an equal superposition of all degenerate states 
except l^i) is created in the lower frame. 

Figure 3] shows the final populations versus the rms 
pulse area A = ttxT for N = 3 degenerate lower states 
for couphngs chosen to satisfy Eqs. ((30)) (upper frame) 
and ([55)1 (lower frame). As follows from Table |TT1 an 
equal superposition of all degenerate states is created for 
rms pulse area A = IStt; this is indeed seen in the figure 
in the upper frame. For the same value of A in the lower 
frame an equal superposition is created of all degenerate 
states except the initially populated state IV'i)- In both 
frames, there are other values of A for which the same 
superpositions are apparently created; a closer examina- 
tion (not shown) reveals that for these other values of the 
rms pulse area the created superposition has almost, but 
not exactly, equal components. 

Figure [5] displays the time evolution of the populations 
for TV = 3 degenerate lower states and rms pulse area of 
IStt, and two detunings: A = in the upper frame and 
AT — 50.534 in the lower frame. For these pairs of areas 
and detunings. Figs. [3] and [4] have already demonstrated 
that an equal superposition of all degenerate states is 
created. Figure [H] shows that the evolution towards such 
a superposition can be dramatically different on and off 




Detuning (units of l/T) 



FIG. 3: (Color online) Populations vs the detuning Aq for 
A'^ = 3 lower states and xT = 18. The coupling strengths Xn 
are given by Eqs. (I30|l in the upper frame and Eqs. (|33p in 
the lower frame. The system is initially in state {"ipi)- 

resonance. Indeed, for A = (upper frame) the nonde- 
generate upper state receives considerable transient pop- 
ulation, which would lead to significant losses if this state 
can decay on the time scale of the pulsed interaction. In 
strong contrast, off resonance this undesired population 
is greatly reduced (lower frame), and still the desired 
equal superposition of the degenerate states emerges in 
the end. We have verified numerically that for larger de- 
tunings this transient population continues to decrease, 
e.g. for AT = 142.198 and QT = 30 it is less than 1%. 

To conclude this section we point out that one can 
create any desired superposition, with arbitrary unequal 
populations, in very much the same manner, on or off 
resonance, by appropriately chosing the individual cou- 
plings, while still maintaning particular values of the 
overall rms pulse area. Tuning on resonance gives the 
advantage of smaller pulse area required, whereas tuning 
off resonance (with larger pulse area) provides the advan- 
tage of greatly reducing the transient population of the 
possibly lossy common upper state. 

VI. MULTISTATE LANDAU-ZENER MODEL 

As seen in Table U the propagator parameter a for the 
LZ model, a = exp (— 7rx^/4C) , cannot be equal to 
or 1 or —1, but may approach or 1 arbitrarily closely. 
However, it is always positive and cannot approach the 
value of —1; hence the LZ model is unsuitable for unitary 
operations within the degenerate manifold, in contrast to 
the resonance and RZ models discussed above. Still, the 
present multistate LZ solution represents an interesting 
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FIG. 4: (Color online) Final populations versus the rms pulse 
area x for A'^ = 3 degenerate lower states and detuning AT = 
50.534. The coupling strengths Xn are given by Eqs. (|3Up 
in the upper frame and Eqs. (I33|) in the lower frame. The 
system is initially in state 
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FIG. 5: (Color online) Populations versus time for A'' = 3 
lower states and rms Rabi frequency xj^ = 18. The coupling 
strengths x-n, are given by Eqs. (|30|l . The detuning is A = 
in the upper frame and AT = 50.534 in the lower frame. The 
system is initially in state IV'i)- 



and important addition to the available LZ solutions (see 
and references therein). 



states It/'at+i) and j^/^n) at the crossing of their energies. 



The Demkov-Osherov model 



2. The degenerate case 



The present multistate LZ model complements the 
Demkov-Osherov (DO) model [1^, wherein a slanted en- 
ergy crosses TV parallel nondegenerate energies. In the 
DO model, the exact probabilities Pn^m have the same 
form — products of LZ probabilities for transition or no- 
transition applied at the relevant crossings — as what 
would be obtained by naive multiplication of LZ prob- 
abilities while moving across the grid of crossings from 
IV'n) to \'ipm)i without accounting for phases and interfer- 
ences. For example, if the states ji/'n) [n — 1,2, . . . , N) 
are labeled such that their energies increase with the in- 
dex n, and if the slope of the slanted energy of state 
l^jv-i-i) is positive, the transition probabilities in the DO 
model are 



p 

^ n — >m 
P 

P 

Pn+i^n+1 



PnQn+lln+2 ' ' ' <Zm- 

(n > m), 

Qn, 

Pnln+lln+l ' ' ' QN , 
qiq2 ■ ■ -qn-lPn, 

qiq2 ---qN, 



(46b) 
(46c) 
(46d) 
(46e) 
(46f) 



where (/„ = cxp (— 7rx^/2C) is the no-transition probabil- 
ity and pn — 1 — q-n is the transition probability between 



The present multistate LZ solution provides the tran- 
sition probabilities for the special case when all parallel 
energies are degenerate, which cannot be obtained from 
the DO model. 

a. The propagator The elements of the transition 
matrix for our {N + l)-state degenerate LZ problem are 
readily found from Eq. (PO)) to be 



7-7- XnXm 

^m,n — 2 V 



U - 1 - AH 
X 



^) (to, n — 1, 



, N; m 7^ n) , 
(47a) 



1-e"'^) {n^l,...,N) 



Un,N+i = —b {n = 1, 
X 

UN+l,n = — — b* 

X 



{n = l 
Un+i,n+i = e 



(47b) 
(47c) 
(47d) 
(47e) 



with A = 7rxV4C and \b\^ = 1 - e^^A. 

b. System initially in the nondegenerate state When 
the system begins initially in the nondegenerate state 
|?/'jv+i), with the tilted energy, the system ends in a co- 
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herent superposition of all states with populations 



P„ = :i|(l-e-2A) („ = i,...,7V), (48a) 



N+l 



-2A 



(48b) 



In the adiabatic limit A 1 the population is distributed 
among the degenerate states according to their couplings, 
whereas the initially populated state \^pN+l) is almost 
depleted, Pn+i ~ 0. For equal couplings, all degenerate- 
state populations will be equal, P„ w 1/7V. In the op- 
posite, diabatic limit A ^ 1 the population remains in 
state \ipN+i) with almost no population in the degenerate 
states. 

c. System initially in a degenerate state When the 
system is initially in an arbitrary degenerate state 
at the end of the evolution the populations are 



Xt 



(1- 



-,-A^ 



(49a) 



Pn - ^(1-e-^) (n = l,...,iV;n^09b) 



N+l 



= ^1(1- 



X 



(49c) 



In the adiabatic limit A 3> 1 and for equal couplings, the 
populations will be 



Pr 



1 

N 



(50a) 



Pn « ^ in = l,...,N;n^i), (50b) 



1 

N' 



(50c) 



Obviously, Eqs. (|48|) and ((49|) cannot be reduced to the 
DO solution (I46p . which implies that the non-degeneracy 
assumption in the DO model is essential. 

Figure [6] shows the transition probability for the mul- 
tistate LZ model plotted against the LZ parameter A = 
7rx^/4C. As A increases the populations approach their 
steady adiabatic values ([50]) . Different coherent superpo- 
sitions can be created by choosing appropriate values for 
the couplings Xn- 



allow one to find special values of parameters, termed 
1.0 




2 4 
Landau-Zener Parameter A 



FIG. 6: (Color online) Populations for the degenerate LZ 
model vs the LZ parameter A = 7rx^/4C for N — 3 degen- 
erate states and equal couplings. The system is supposed to 
start in one of the degenerate states \tpi). The arrows on the 
right point the adiabatic values (f50)) . 



generalised tt pulses, for which various types of popula- 
tion transfer can occur, for example, creation of maxi- 
mally coherent superpositions. The RZ solution is par- 
ticularly useful because it allows to prescribe appropri- 
ately detuned pulsed fields for which the dynamics can be 
essentially contained within the degenerate-state space, 
without populating the upper state even transiently, thus 
avoiding possible losses from this state via spontaneous 
emission, ionization, etc. 

We have analyzed in some detail also the multistate 
Landau-Zener model, which complements the Demkov- 
Osherov model in the case of degenerate energies. 

The presented analytical solutions and general proper- 
ties have a significant potential for manipulation of mul- 
tistate quantum bits in quantum information processing, 
for example, in designing arbitrary unitary gates. 



VII. CONCLUSIONS 

In this paper we have described a procedure for deriv- 
ing analytical solutions for a multistate system composed 
of N degenerate lower states coupled via a nondegenerate 
upper state with pulsed interactions of the same tempo- 
ral dependence but possibly with different peak ampli- 
tudes. The multistate resonance and Rosen-Zener solu- 
tions have been discussed in some detail because they 
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